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Correlations between intrinsic ellipticities of galaxies are a potentially important systematic error 
when constraining dark energy properties from weak gravitational lensing (cosmic shear) surveys. 
In the absence of perfectly known galaxy redshifts some modeling of the galaxy intrinsic alignments 
is likely to be required to extract the lensing signal to sufficient accuracy. We present a new model 
based on the placement of galaxies into dark matter halos. The central galaxy ellipticity follows the 
large scale potential and, in the simplest case, the satellite galaxies point at the halo center. The 
two-halo term is then dominated by the linear alignment model and the one-halo term provides a 
motivated extension of intrinsic alignment models to small scales. We provide fitting formulae for 
the spatial projected source power spectra for both intrinsic-intrinsic (II) and shear-intrinsic (GI) 
correlations. We illustrate the potential impact of ignoring intrinsic alignments on cosmological 
parameter constraints from non-tomographic surveys, finding that as could be underestimated by 
up to the size of the current l-cr error bar from cosmic shear if very small scales are included in the 
analysis. Finally, we highlight areas of interest for numerical simulations of dark matter clustering 
and galaxy formation that can further constrain the intrinsic alignment signal. 

PACS numbers: 98.80.Es, 98.62.Gq, 98.62.Sb 



I. INTRODUCTION 

Gravitational lensing causes the images of distant 
galaxies to be distorted as the light travels through space- 
time that has been bent by intervening matter. Two 
galaxies that are physically close will therefore appear 
preferentially aligned with each other because their light 
travels through similar regions of curved space-time, so 
their images are both distorted in a similar direction. 
This "cosmic shear" seems to be one of the most promis- 
ing probes of the nature of the mysterious dark en- 
ergy [l|, [2| that apparently dominates the energy budget 
of our universe. 

However, the lensing distortions are extremely small 
compared to the intrinsic ellipticities of galaxies, and can 
only be detected by averaging the ellipticities of many 
spatially localized galaxies. If galaxies are randomly ori- 
ented in space then any intrinsic ellipticities should aver- 
age to zero allowing a detection of the gravitational lens- 
ing signal. But because we believe galaxies form inside 
large-scale gravitational potentials, coherent tidal effects 
might be expected to align the intrinsic galaxy elliptici- 
ties. 

There are two physical mechanisms that have been 
widely considered for inducing alignments of galaxy ori- 
entations. First, there may be alignment of angular mo- 
mentum vectors of galaxies that formed in the same ini- 
tial tidal field due to tidal torquing [see e.g. Q. Sec- 
ond, there may be a coherent tidal stretching of galaxy 
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shapes due to the large scale gravitational potential |5| 
or anisotropic accretion along filaments Q. The align- 
ment of spin axes is believed to dominate for rotationally 
supported spiral galaxies while the tidal stretching is be- 
lieved to be the dominant effect for elliptical galaxies. 

Using linear tidal torque theory, Rcfs. 0, [|[ showed 
that the correlation function for the galaxy ellipticity 
orientations is nonzero only at second order in the tidal 
tensor and is therefore small except at small galaxy sepa- 
rations. The tidal stretching of elliptical galaxies can be 
quantified analytically by assuming that the intrinsic el- 
lipticity of galaxies is proportional to the curvature of the 
primordial large scale potential [aL This is often referred 
to as the linear alignment model [9j. Since the ellipticity 
correlation function of tidally stretched elliptical galax- 
ies is linear in the tidal tensor it could potentially have 
a significant amplitude over many tens of megaparsecs. 

The size of these effects has been estimated using 
numerical simulations, in which galaxies are effectively 
pasted into the simulation with directions corresponding 
to the orientation of individual dark matter halos, or the 
angular momentum vector [h], EH, Ell EH • This has con- 
firmed the result that the most significant correlations 
come from alignments of galaxy light with the orienta- 
tion of the dark matter halo [13| . 

On slightly smaller scales some common trends have 
been observed, both in simulations and observations, re- 
garding the alignment of sub-halos or satellite galaxies 
within their parent halo or cluster/group. First, satel- 
lites tend to be aligned with the radius vector of their 
host halo [3, El, El El El, El, HJ- Second, the satel- 
lites are found to be preferentially located near the major 
axis of the halo [U, lU, EI HE H] ■ Third, the satellites 
are aligned with the central galaxy [ll| [24[ ■ See Ref. [25[ 
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for the most recent observations of these effects. 

The alignment of physically neighboring galaxies is ex- 
pected to produce a spurious increase in the cosmic shear 
galaxy alignment signal. This is often referred to as the 
II (intrinsic-intrinsic) alignment. Ref. 0] pointed out 
that an additional contamination can occur due to the 
simultaneous alignment of a galaxy with a nearby mat- 
ter distribution, and the gravitational lensing of a more 
distant galaxy by the same matter distribution. This 
gravitational-intrinsic (GI) alignment produces a spuri- 
ous decrease in the cosmic shear signal, since the two 
galaxies now point in opposite directions. The closer 
galaxy points towards the mass clump, whereas the 
more distant galaxy is stretched tangcntially around the 
clump. A small-scale contribution to the GI effect even 
comes from the gravitational lensing of distant galaxies 
by a galaxy's own elliptical halo (26J. 

Observations of the II effect have been carried out us- 
ing low redshift galaxy samples by SuperCOSMOS [27j . 
COMBO-17 [H[ and SDSS (H,^. The raw GI effect 
is harder to measure, due to the large range of redshifts 
involved, and potential contamination by cosmic shear it- 
self. However the alignment between galaxies and neigh- 
boring mass has been measured by using the distribution 
of galaxies themselves as a tracer of the mass distribu- 
tion mm 

The methods for removing both II and GI intrinsic 
alignment (IA) effects as a systematic error in cosmic 
shear measurements can be classified in two categories of 
"nulling" or "modeling" [32j . "Nulling" methods down- 
weight selected parts of the data in such a way as to 
remove the intrinsic alignment signal. For II correlations 
this can be done by removing galaxy pairs that are close 
in redshift (and on the sky) [281 [33LI33. [35|. while the GI 
signal can in principle be removed with a particular linear 
combination of tomographic shear power spectra [3(| . 

The "modeling" technique instead specifies a param- 
eterized model of the intrinsic alignments and then 
marginalizes over the instrinsic alignment parameters 
when inferring cosmological parameter constraints (i.e. 
the systematic error is subtracted and reduced to a sta- 
tistical error) [33|, [37j • This has the advantage of not re- 
quiring precise galaxy redshift measurements, but has the 
disadvantage of being potentially sensitive to the choice 
of intrinsic alignment model and the number of parame- 
ters used in the model (with more parameters potentially 
leading to larger degradations in the inferred cosmolog- 
ical constraints [13|)- Recent work has also shown that 
it may be possible to "self-calibrate" the IA signal using 
the cross-correlation between the ellipticity and galaxy 
density fields in the same survey [38|. 

Ideally a sufficiently good prediction of the IA signal 
could be made from simulations and subtracted from the 
data to leave the pure lensing signal. However current 
simulations have some way to go before they will make 
reliable predictions, due to the finite resolution of a large 
box size, and the difficulty of including baryonic material. 
And in any case it will be necessary to have predictions 



for a range of cosmological models and thus an analytic 
fitting formula would be very helpful. The linear align- 
ment model is likely to hold on large scales, whereas it 
takes no account of non-linear growth of structure. An 
ad-hoc attempt to rectify this was indicated in Ref. [3l| , 
and used in Ref. (37j . in which the linear theory mat- 
ter power spectrum in the linear alignment model was 
replaced by the non-linear matter power spectrum. We 
refer to this as the nonlinear alignment model (NLA) 
here. In this paper we construct an improved model for 
the IA power spectra based on the halo model. 

The halo model of galaxy clustering [39|, H(j has proved 
a surprisingly successful predictor of galaxy clustering 
statistics [see Hi], for a review] . In this picture the uni- 
verse consists of dark matter halos that are clustered ac- 
cording to linear theory (in the case of two-point function 
predictions). Each halo has a mass drawn from a mass 
function, and a density profile. These are usually taken 
from average properties of n-body simulations. Galaxy 
positions are drawn from the resulting dark matter distri- 
bution. The two-point correlation function of galaxy po- 
sitions then constitutes a "two-halo" term, arising from 
the correlations between the positions of two different ha- 
los, and a "one-halo" term arising from pairs of galaxies 
that reside in the same halo. 

This paper is organized as follows. We describe our 
model for central and satellite galaxy ellipticities in Sec- 
tion |TTJ In Section IIIII we derive the 3D auto power 
spectra of the projected intrinsic galaxy ellipticity dis- 
tribution and the cross-power spectra with the matter 
distribution. We derive the intrinsic alignment contri- 
bution to the angular shear power spectra in Section HVl 
and show how constraints on og could be biased if the in- 
trinisc alignment correlations were modeled incorrectly. 
In Section [V] we draw conclusions about this new satel- 
lite contribution to the intrinsic alignment signal and de- 
scribe some questions for future numerical simulations 
to address in order to refine our model. We describe the 
normalization of the 3D ellipticity power spectra account- 
ing for spatial clustering of the galaxies in Appendix [A] 
and give details on our analytic model calculation in Ap- 
pendix [B] 

Unless otherwise stated, we assume a fiducial cosmo- 
logical model with a s = 0.8, O m = 0.3, flk = 0, n s = I, 
and H = 71.9 km s" 1 Mpc" 1 . 



II. MODEL FOR GALAXY ELLIPTICITIES 

We are concerned with modeling the correlations in 
the intrinsic projected orientations of galaxies that could 
mimic those induced by gravitational lensing. There are 
three effects to consider when constructing such a model: 
the 3-D shapes of galaxies, projection effects, and the 
relation of the projected ellipticity to the shear defined 
in gravitational lensing studies. 

Elliptical galaxies can generally be described by 3-D 
ellipsoids, which are further found to be mostly pro- 
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FIG. 1: Cartoon depiction of the galaxy alignments within 
and between halos. The intrinsic alignment correlations 
within a single halo come from the radial alignment of the 
satellites (green lines) . The correlations between separate ha- 
los are dominated by the correlation of the central galaxies in 
each halo, depicted as red lines. 

late [42j]. Assuming the two minor axes of the ellipsoid 
have equal lengths on average, we can then approximate 
elliptical galaxies as sticks corresponding to the length 
and orientation of the major axis. In this paper we as- 
sume that all halo satellite galaxies can be described in 
this way. If the stick makes an angle 9 with respect to the 
line-of-sight we relate the observable projected intrinsic 
cllipticity of the galaxy to the length of the stick 7 as, 

17^^-1 = 7 sinfl. (1) 
a + 

where a and b are the major and minor axes of the 2-D 
projected ellipticity of the galaxy. We follow the usual 
convention and define the complex projected cllipticity 

7 / (r) = |7 / (r,0)|e 2 ^^7i(r)+*72(r). (2) 

Figure[T]is a schematic picture of our stick model show- 
ing how we populate spherical dark matter halos with 
central and satellite galaxies. We describe the separate 
intrinsic alignment models for centrals and satellites in 
the next two sub-sections. 



A. Satellite galaxies 

We build a model for intrinsic alignments of satellite 
galaxies by populating a spherical dark matter halo with 
galaxies, with a number density that follow the density 
profile of the halo. Our "basic" model for satellites as- 
sumes that all satellite galaxy ellipticities are aligned 
with the radius vector the halo, in three dimensions. This 
model is inspired by the idea that tidal forces within the 
halo are responsible for the intrinsic alignments (20| . 

In a flat-sky coordinate system with the origin in the 
center of the halo, the z axis aligned with the line-of-sight 



direction, and polar angle 9, the 3-D density-weighted 
projected ellipticity of galaxies in the halo is (using the 
same notation as Ref. Q) 

7 / (r, ra, c) = 7 / (r,m, c)N g u(r\m, c) (3) 
= 7(r,m,c)e 210 sin 9 N g u{v\m,c) (4) 

where 7(7*, to, c) is the magnitude of the projected elliptic- 
ity at radius r in a halo of mass to, and c, which denotes 
the concentration of the NFW halo profile Pnfw [HI- 
We follow Ref. 0] in working with the density-weighted 
cllipticity. In general, j(r, to, c) may be an arbitrary 
(positive-valued) function of the radial position within 
and mass of the halo, and may vary from halo to halo. 
However, in all the numerical calculations we perform be- 
low we set 7(7", to, c) = 0.2, independent of position, mass 
and concentration. 

An illustration of the cllipticity correlations induced 
by this radial stick model inside a single halo is shown 
in Fig. O An example galaxy is shown in green in the 
middle right of the figure, and the ellipticities of three 
concentric circles of galaxies are shown around this using 
black lines. Note that these lines vary in length due to 
the density weighting and there are no galaxies belonging 
to this halo outside the virial radius, denoted by the black 
circle. Components of the black lines parallel and per- 
pendicular to the green line are shown in red and blue 
respectively. The cartoon shows that this radial align- 
ment model gives a characteristic shape to the one-halo 
correlation function (as usual £±(#) = £++(#) ± £xx(#) 
where £+,x(#) = {{7+,x (d')l+,x + 0))e')directions of e, 
where 7+ and 7 X are the components of 71 and 72 in the 
coordinate system aligned with 6 - see also Eqn. (IA1[) ). 
The red lines clearly dominate for the central circle, caus- 
ing a positive correlation. The next circle is dominated 
by perpendicular components shown in blue, producing 
a net negative correlation. Finally the largest circle con- 
taining any galaxies within the virial radius will produce 
a positive correlation. See the black line in Fig. for an 
example correlation function (averaged over all positions 
of the green line in Fig. [2] as usual). We computed this 
correlation function both using the spherical harmonic 
approximations used throughout this paper and detailed 
in the Appendices, and also using a Monte Carlo sim- 
ulation of galaxies placed into an NFW halo. We find 
excellent agreement. The one- halo satellite-satellite term 
in the left panel of Fig. 0] shows the E-mode power spec- 
trum constructed from the Fourier transforms of £± in 
Fig. GH (according to, e.g., Eqn. (|A1[) ). 



B. Central galaxies 

We use the linear alignment model [1, Qj| 5H to spec- 
ify the intrinsic ellipticities of central galaxies in our halo 
model. Central galaxies are defined such that each dark 
matter halo has one central galaxy that resides exactly 
at the center of the spherical halo mass distribution. We 
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FIG. 2: Depiction of intrinsic alignment correlations of satel- 
lite galaxies inside a single halo of mass 10 14 M Q /h. The green 
line is ellipticity of a typical galaxy that we may be correlat- 
ing with all the other galaxies in the cluster. The black lines 
show density weighted ellipticities of other cluster galaxies at 
3 fixed radii from the green galaxy. The red and blue lines 
are the components of the black lines denoting the degree of 
correlation with the green galaxy. Red indicates positive cor- 
relation while blue indicates negative correlation. The total 
correlation from the 3 circles can be obtained by adding the 
red and blue lines by eye - which makes it very clear that the 
smaller circle is strong positive correlation while the middle 
circle has a weak negative correlation and the largest circle 
again has a small positive correlation. 




6 [arcmin.] 

FIG. 3: Effect of a distribution in radial alignment angles 
of the satellite galaxy ellipticities on the one-halo correlation 
function. The upper lines (black/solid and blue/dashed) are 
the two ellipticity correlation functions for a single halo of 
mass M = 10 14 h^ 1 Mq and concentration = 5 at z = 0.3 
with perfect radial alignments of the satellites. (Absolute 
values are shown.) In yellow (gray) are the same correlation 
functions when the radial alignments angles of the satellites 
are perturbed by the distribution from Ref. [l!| (also given in 
Eqns. [16] and I17|l . The dotted lines show the black and blue 
upper lines rescaled by the multiplicative factor 7s Ca i G = 0.21 2 . 
The dotted lines almost perfectly overlay the correlation func- 
tions with the distribution in alignment angles. 



III. INTRINSIC ELLIPTICITY POWER 
SPECTRA 



assume that central galaxies have the same ellipticity ori- 
entation as their parent halo and that halo ellipticities 
are determined by the large-scale density perturbations 
so that [eq. 13 in Ref. 

jl r = -^( V ^ V ?' 2V ^) S ^ ^ 

where iSf^p] is the primordial potential with small-scale 
perturbations smoothed out and C\ is a normalization 
constant. Unless otherwise stated we take C\ — 5 x 
10 14 (h 2 M Q Mpc 3 ) -1 . This is motivated by compari- 
son with Fig. 2 in Ref. 0] who normalise to SuperCOS- 
MOS Si and that fits well with the SDSS L4 points 
of Ref. [29j as illustrated in Ref. pT^ . 

We can recover the linear alignment model in our 
framework by setting the Fourier transform of 7 / (r) to, 

C o 

7 J (k) = -j£ sin 2 9 k (cos(2</> fc ), sin(20 fc )) 5 lin (k), (6) 

where <pk is the azimuthal angle of the Fourier wavenum- 
ber k about the line-of-sight, 8k is the polar angle of k 
with respect to the line of sight, Sn n is the linear theory 
density perturbation and the remaining definitions follow 
Ref. [|. 



We construct a continuous intrinsic ellipticity field by 
following the usual halo model assumption (e.g. Ref. |45j ) 
in which we sum the density- weighted ellipticity distribu- 
tion for a single halo over halo masses rrii and positions 
rt, 

7 J (r) = ^-^7 / (r-r l ,TOi)A r 9!l it(r-r 4 |m 4 ) 

n s i 

x ^ 7 f (r-r',m)«(r-r'|m), (7) 

n g 

where u(r\m) = pNFw( r j rn)/m, N g j is the number 
of galaxies in the zth halo, and n g is the mean num- 
ber of galaxies per unit volume (at a given redshift). 
We will assume a deterministic relation between halo 
mass m and concentration c hereafter with c(M, z) = 
{9/(l + z))(M/M*(z))-°- 13 0. 

We then compute the 3-D ellipticity power spectra by 
first Fourier transforming each component of Eqn. ([7]), 

7/(k,m)= /rf 3 r 7 /(r,m)e jk - r (8) 
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where j = 1,2 denotes the shear component. We perform 
this Fourier transform for our model using a multipole 
expansion of the plane waves as shown in Appendix IB 21 
The E and B Fourier modes are defined as 



7i(k)=cos(20 fe )7 1 / (k) 
7 i(k) = sin(2«/> fe ) 7l (k) 



-sin(2^.)7 2 / (k) 
cos(2^) 72 7 (k). 



(9) 



We then define the 3-D ellipticity power spectra, 

<7lT(k)7£(k')) = (2^) 3 ^ 3) (k - k') P EE (k) 
(7^(k) 7 ^(k')) = (2n) 3 S^ (k - k') P BB (k) 



and (Ng(Ng — 1)|tti) are the first and second moments 
of the galaxy number distribution within a halo of mass 
to, and the s and c subscripts denote "satellite" and 
"central" galaxies, respectively. 

The two-halo central-central power spectrum is 
(Eqn. (16) in i), 

Pif c f (fc)=C?iW*)+ 



/,/3l 

X /s(ki)fl in (fci).Fij n (fc 2 ), 



(12) 



<<y*(k)7|,(k')> = (27r) 3 4 3) (k-k')n,^(k). (10) where k 2 = k - ki , C = C lP /D, b g , c is the linear galaxy 

bias for the central galaxy population, and, following 
Ref. 0, f E (k) = (kl - kf) /k 2 . The two-halo central- 
satellite term in the power spectrum is, 



When considered as contributions to the cosmic shear 
signal, P~j X gives the II contribution from the intrinsic- 
cllipticity correlation of two galaxies in the lens plane, 
while P$^i gives the GI contribution from the correlation 
of a galaxy in the lens plane with its surrounding dark 
matter distribution (which is the lens). 

First we consider the case where all the satellite galax- 
ies are pointing directly towards the center of the halo. 
Then we consider the case where there is some random- 
ness in this direction. We then present fitting formulae 
for the above power spectra. 



A. Radially aligned satellites 

Two types of term appear when taking products of 
the form in Eqn. (fTDj) : both galaxies are in the same 
dark matter halo or they are in different halos. These 
are commonly referred to as the "one-halo" and "two- 
halo" contributions to the power spectrum. In addition, 
because we have separate models for central and satellite 
galaxies in halo (and because the probabilities for cen- 
trals and satellites to occupy a halo of mass m may be 
different) we get an additional separation of terms into 
correlations of satellites with satellites, centrals with cen- 
trals, and satellites with centrals. 

By integrating over the joint probablitity distribution 
for 2 halos of masses toi and to 2 to be at positions 17 and 
r 2 with N 9t i and -/V ffj2 galaxies in each halo we arrive at 
the satellite-satellite 3-D II power spectra, 



pEE,lh 

7 7 7 SS 



(JV«(iV'-l)|m) 
amn(m) - — - ^ - 



P?i E 's h (k)= I dmm{mi) 



x 7 2 (ra) \w(k,8k\m)\ 2 

(N> g \mi) 



(N s g \m 2 ) 



x P hh {k\m 1 ,m 2 ), 



x / cZto 2 ti(to 2 ) 



7(toi) |io(A;,0fc|mi)| 



7(m 2 ) \w(k,9 k \m 2 )\ 



(11) 



where P^ L is the halo-halo power spectrum, 
w(k\m) = 7 / (k,TO)/7(m) (sec Eqn. (|B5[) i. (N g \ni) 



pEE,2h 



(k) 



Cip 



D 



PUk) 



dmn(m) * bh(m) 



x 7(771) \w(k,9k\m)\ 



(13) 



where is the halo bias. In principle there could be 
a one-halo central-satellite correlation, but this is zero 
when using the linear alignment model for the central 
galaxies and assuming that the satellite galaxy elliptici- 
ties are uncorrelated with the central galaxy ellipticity. 

Because we have normalized by the comoving den- 
sity of galaxies n g (z), these power spectra are only very 
weakly dependent on the model rcdshift distribution. 
We have verified that the weight in the II integrand, 
mn{m, z) (N g \m) jn g {z) has an r.m.s. error of less than 
10% between redshift distributions with median z ~ 0.6 
and ~ 1.7. Therefore we consider only a single rcdshift 
distribution for the remainder of this paper (which is the 
model from Ref. (4(| with limiting magnitude in R of 
26.). Changing the relative contribution of spiral and el- 
liptical morphologies to the satellite population, and as- 
suming that spirals have zero intrinsic alignments, does 
affect the power spectrum by diluting the IA signal. Us- 
ing ad-hoc models for the halo mass-dependent ratio of 
spiral and elliptical satellites we find the dominate effect 
is to simply rescale the amplitude of the power spectra. 

We show in Appendix IB 21 that 7(3 s = 0; implying 



pBB.lh _ pBB.Ih 



pB ( ,2 _ q rphjg f rj ows directly 

7 > cs 



from the simple, and separable, ^-dependence in 7 7 (r). 
We can understand this more directly by noting that all 
the satellite galaxy ellipticities are perpendicular to the 
boundary of the halo (assumed to be truncated at the 
virial radius). From, e.g., Ref. we know that this 
boundary condition is incompatible with a nonzero B- 
mode 1 . The linear alignment model B-mode, P? T B ' 2h was 



1 Although, this condition does not distinguish "pure" E modes 
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FIG. 4: All nonzero contributions to the II and GI E-mode power spectra in the spherical, radially aligned satellite, halo model. 
For comparison we plot the "nonlinear alignment model" in which the non-linear matter power spectrum is used within the 
linear alignment model. 



derived in Ref . Q , but is much smaller than the E-mode 
and we do not consider it further. 

Some of the two-halo terms we have to consider for an 
ensemble of halos are demonstrated in cartoon form in 
Fig. Q] The dominant contribution to the two-halo term 
is the correlation of the central galaxies. This makes 
sense intuitively if we think of averaging the satellite el- 
lipticity correlations by eye in Fig. [TJ 

The nonzero GI power spectra are 



P &,s = I dmn(m) 



m [N g \m) 



P n g 

x 7(771) \w(k,9k\m)\ u{k\m) 



7(toi) \w(k,9 k \mi)\ 



x / drri2 n{m2) — u(k\m2) 
J P 

x P hh (k\mi,m 2 ) 



p2h 



D 



flin(fc), 



(14) 



where P$~i c is the same as Eqn. (18) in Ref. @. We as- 
sume that the halo-halo power spectrum P^h is given by 
the linear theory power spectrum Py ln multiplied by the 
halo bias for each mass. Again, there could be a one- halo 
central galaxy term, but it is zero for the assumptions we 
have made in this halo model. 

We show the various one-halo and two-halo compo- 
nents of the E-mode II and GI power spectra in Fig. |U 
The total power spectra are well-described by only con- 
sidering the satellite contribution to the onc-halo term 



from ambiguous modes as defined in Ref. [47 
tion is not important for this simple model. 



But, this distinc- 



and the central contribution to the two-halo term. We 
will neglect the other cross-terms shown in Fig. 0] in our 
remaining analysis. We see that the one-halo contribu- 
tions to both the II power spectrum becomes important 
around k = lh Mpc -1 . This fits quite well with the in- 
crease in power in the non-linear matter power spectrum. 
The amplitude on small scales is similar to the nonlinear 
alignment model for the II power spectrum, and signifi- 
cantly larger for the GI power spectrum. 



B. Distribution in radial alignment angles 

In this section we consider the effect of a distribution 
in the angle between the galaxy major axis and the radial 
vector for satellites in a halo. If f is a unit radial vector 
in the halo and e is a unit vector denoting the orientation 
of the galaxy major axis, then we consider a distribution 
in the angle (3 where, 



cos /3 



(15) 



By default, we adopt the distribution found in Ref. [l| 
(their Eqn. (5)) from N-body simulations, 



P(f3) = sin(/3) 



.4 



B + A/5 



cos 4 (3 + B 



(16) 



with A = 2.6 and B = 0.6. Once we allow a misalignment 
between the galaxy major axis and the halo radial vec- 
tor then there is a second angle to consider that denotes 
rotations about the axis defined by f, which we label rj. 
We assume a uniform distribution, 



(17) 



We derive the satellite ellipticity as a function of (3 and 
rj in Section UTTl We have also performed a Monte Carlo 
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C. Fitting functions 

We show the II and GI satellite power spectra at sev- 
eral redshifts in Fig. [5] with 7 sca i e set to unity. The ampli- 
tudes of both spectra decrease with increasing redshift, 
which is largely a consequence of the decreasing fraction 
of satellite versus central galaxies as the redshift increases 
(so the one- halo signal is essentially diluted). The peaks 
of the spectra also shift to smaller scales with increasing 
redshift as the objects move further away. 

We find the one-halo II and GI power spectra at a 
given redshift can be reasonably fit with the following 
3-parameter functions, 



FIG. 5: Ellipticity suppression factor -y sca i for different levels 
of random perturbations /3 in the radial satellite alignments. 
7scaio is shown as a function of the Gaussian width of P(cos /3) 
with the result from Ref. nyl overlaid. 
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simulation of a single halo to more robustly study the 
effect of this distribution in alignment angles on the one- 
halo correlation functions. 

The resulting correlation function is shown by the 
lower pair of lines in Fig. We find the main effect 
of adding a distribution in radial alignments is simply to 
reduce the amplitude of the correlation functions with 
respect to the case with perfect radial alignemnts. For 
the distribution in Eqn. (|16p , we get a multiplicative am- 
plitude reduction of 



We show the evolution of pi as a function of redshift 
in Fig. We fit each p t (for i = 1, 2, 3) for both II and 
GI fitting functions with the model, 

Pi(z) = qii exp (q l2 z qa ) . (20) 

The best fit values for the qij are given in Table Q] for the 
II spectra and in table [TT1 for the GI spectra (note that qn 
has units of /iMpc -1 ). The fit functions are compared 
to the input spectra in fig. [7] 



%'caic = 0-21 2 (18) 

independent of halo mass. The dot-dashed lines in 
Fig. ([3]) show the correlation function without a distribu- 
tion in alignment angles, scaled down by this factor. It is 
in good agreement with the full calculation. Unsurpris- 
ingly this suppression factor can also be obtained more 
quickly by integrating over (3 and v for a satellite at a 
given polar angle 9 in the halo. 

We have also performed simulations using a Gaussian 
distribution in the radial alignment angle (3 with zero 
mean and a varying width a p. The amplitude reduction 
as a function of ap is shown in Fig. ([5]) where we find 
that 073 ~ 0.7 gives similar correlation functions to those 
obtained using the distribution in Eqn. ([TBI . 

When P(/3) has finite width, a nonzero B-mode can 
be generated for 7J. We have computed this amplitude 
via numerical integration of fg in Eqn. (|B8[) (and us- 
ing Eqn. IB4[) with a Gaussian distribution for P(j3) and 
find the B-modc amplitude is several orders of magnitude 
smaller than the E-mode amplitude independent of the 
width of the P{0) distribution. We therefore ignore the 
B-mode power spectra in the following results. 



TABLE I: Fitting function parameters for pi in the II power 
spectrum fit. 



param. index 


qn 


qa 


qn 


1 


0.09939 


3.718 


0.3475 


2 


1.931 


1.061 


0.7484 


3 


6.082 


0.1045 


0.613 


TABLE II: Fitting 


function paramet 


ers for pi in 


the GI power 


spectrum fit. 








param. index 


qn 


q%i 


qn 


1 


0.01867 


6.924 


0.3725 


2 


1.989 


1.081 


0.6816 


3 


4.232 


0.1748 


0.481 



We then have 9 parameters each to describe the red- 
shift dependent II and GI one-halo power spectra, since 
qn is degenerate with 7 sca i e (as well as the fraction of 
spiral to elliptical galaxies in a halo as mentioned in Sec- 
tion lllip . Combined with the free amplitudes for the two- 
halo terms, we have 20 parameters in total to describe 
the intrinsic alignment power spectra. 
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FIG. 6: Comparison of the II and GI one-halo power spectra and their fitting functions at 9 different redshifts. The dashed 
(red) lines show the best fit of Eqn. {19} , while the dot-dashed (green) lines show the best fits of the z-dependent function in 
Eqn. flUJ. 




IV. COSMOLOGICAL IMPLICATIONS 

We use our model to compute all contributions to a 
set of angular cosmic shear power spectra with redshift 
distribution 



signal is given by the halo model with fiducial parameter 
7scalc = 0.21. 



A. Angular shear power spectra 



dN 
dzdfl 



exp 



(21) 



We consider the amplitude of the two-halo term, and 
the amplitude of the one-halo term as parameterized by 
the ellipticity suppression factor, 7 sca io (Eqn. (fT8|) ). that 
we use a proxy for the degree of radial alignment of the 
satellites in a halo. By default we assume the "true" IA 



For comparison with weak lensing power spectra, we 
need to integrate the 3-D projected ellipticity model over 
the lensing source distribution 



dxf a (x)r(xo,x), 



(22) 



where f a (x) is the number of sources in sample a per 
unit comoving distance normalised so J f a (x)dx = 1- 
Therefore, assuming Limber's approximation, the shear 
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FIG. 8: Angular shear power spectra for a fiducial deep survey with 10 tomographic redshift bins and a Gaussian photometric 
redshift uncertainty 0.05(1 + z). The solid line shows the usual power spectra from lensing alone (everywhere positive). The 
dashed lines show the II contribution (everywhere positive). The dotted lines show the GI effect (everywhere negative) and 
the dot-dashed lines show the total (sum), which can in general oscillate positive and negative depending on the balance of 
competition between the three effects. 



angular power spectra can be written 

Q(a/3) = C?{af3) + ^l(a/3) + ( 2 ^) 

where the terms are sourced by the lensing distortions 
due the the matter power spectrum and the intrinsic 
alignment terms calculated in the previous sections 

Jo X 

Jo X 

ClU = r^#^^(fc;x)dx (24) 
Jo X 

where k = i/x- 

We first consider a relatively deep survey with a = 
2, b = 1.5, zq = 0.64 with 10 tomographic redshift bins 
and a photometric redshift uncertainty of 0.05(1 + z). 



This is more relevant for a surveys such as LSST 2 , Eu- 
clid 3 or JDEM 4 . A subset of the resulting angular power 
spectra are shown in Fig. [5] The one-halo contribution 
is clearly seen as the sharp rise in power in II and GI at 
small scales. For this photometric redshift uncertainty 
the II term is quite comparable or subdominant to the 
GI term, even in the autocorrelation bins. The effect of 
intrinsic alignments is clearly largest at very low redshift, 
as expected. There is a competition between the GI, II 
and GG terms that causes the total power spectra for 
the lowest bins to oscillate above and below zero (only 
the absolute magnitudes are shown, but the sign can be 
identified by recalling that the II and GG power spectra 
are positive and the GI power spectrum negative). 

Fig.[5]shows the angular power spectra for a CFTHLS- 



2 www.lsst.org 

3 dune-mission.net 

4 jdem.gsfc.nasa.gov 
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FIG. 9: Angular shear power spectra showing contributions from GG, GI, and II terms for a range of different intrinsic 
alignment models for a fiducial medium-deep survey. The dashed line shows the power spectra for the linear alignment model, 
or equivalently the two-halo term. The dotted line shows the one-halo term only. In light/green dashed is the linear alignment 
model using the non-linear matter power spectrum in place of the linear matter power spectrum. The solid line is the model 
proposed in this paper, consisting of the one and two halo terms. The power spectra in the absence of intrinsic alignments is 
shown by the dot-dashed line. The bottom right panel shows the difference between the total angular power spectra for the 
various models, as compared to the no intrinsic alignments case. 



like survey 5 with a = 0.836, b = 3.425, z = 1.171 in 
Eg. |2"T1 (taken from Ref. [ill). For comparison with exist- 
ing results we use a single redshift bin and a photometric 
redshift uncertainty of 0.1(1 + z). The signal from lens- 
ing alone is shown in the top left panel. The top center 
panel shows Cf 1 for various different models for intrin- 
sic alignments. We see that the one-halo term becomes 
important at multipolc numbers greater than I ~ 500 
and matches surprisingly well to the NLA result. Includ- 
ing either the one-halo term or the NLA boosts the GI 
contamination at I ~ 3000 by an order of magnitude rel- 
ative to that from the linear-alignment model alone. The 
II contribution is shown in the top center panel, where 
we see that the one-halo term dominates on all scales 
where the intrinsic alignment contribution is significant. 
At the smallest scales considered the one-halo term is a 
few times larger than the NLA model and an order of 
magnitude larger than the linear alignment model. How- 
ever both are much smaller than the GI contribution at 
small scales. 

The total effect is shown in the lower panels of Fig. [9] 
There is an overall suppression of the power for all the 
intrinsic alignment models considered, relative to the 
alignment-free power spectrum. This is due to the wide 
range of redshifts covered by the survey we consider. For 
a tomographic survey different contributions would be 
more important for different pairs of tomographic bins, 
as illustrated for the NLA in Ref. [37|- Overall our halo 
model result is qualitatively similar to that from the 



5 www.cfht.hawaii.edu/Scicnce/CFHLS/ 



NLA. The suppression is roughly a factor of two larger 
at the smallest scales considered. On intermediate scales 
the GI suppression is outweighed by the boost in power 
from II, which is more dominant at larger scales than GI. 
The total power spectrum including only the one-halo in- 
trinsic alignment term is actually fractionally larger than 
the no IA power spectrum at intermediate scales I ~ 200. 
Then at larger multipole numbers I > 500 the GI term 
dominates again causing a net reduction in power. 



B. Bias on erg 

Cosmological parameter constraints from cosmic shear 
can have systematic biases if our model for the intrin- 
sic alignments contribution to the shear power spectra 
is wrong or insufficiently flexible. In this section we il- 
lustrate the potential biases on a a for different assumed 
intrinsic alignment parameter values. To find the bias, 
we calculate the fiducial angular power spectra, the fitted 
angular power spectra and the Fisher matrix for the fidu- 
cial spectra for the assumed survey. These are combined 
using Eqn. (21) of Ref. to find the bias on cosmo- 
logical parameters. In this case we are varying only one 
cosmological parameter, erg- We use a maximum I value 
of 2e4, 12 galaxies per square arcminute and 100 square 
degree survey area. 

In the left hand panel of Fig. ([TO)) we plot the bias 
in the inferred <j% value as a function of the ampli- 
tude of the linear alignment model, or two-halo term. 
The fiducial model has our default two-halo term am- 
plitude of Ci = 5 x 10" 14 (/i 2 M Mpc 3 ) _1 and uses 
the fit functions for q given in the previous section, with 
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FIG. 10: Matter power spectrum normalisation parameter erg values that would be inferred when assuming the wrong model 
for the intrinsic alignment power spectra, as a function of C\ (left) and 7scaic (right). The intersection of the black vertical and 
horizontal lines shows the fiducial model with erg = 0.8, Ci = 5 x 10~ 14 (h 2 Mg Mpc 3 ) and 7 sca i = 0.21. The solid (red) 
line shows the inferred erg values using this model with different assumed C\ or 7 sca i c . The dashed line in the left-hand panel 
shows the result when the non-linear alignment model is fitted to the simulated halo model power spectra. The dashed line 
in the right-hand panel shows the result of fitting the halo model to a nonlinear alignment model set of simulated data. The 
dot-dashed line in the left-hand panel shows the result of fitting the linear alignment model to simulated halo model data. The 
cross denotes the erg value that would be obtained on ignoring intrinsic alignments in a naiive fit to data including intrinsic 
alignments from this halo model (i.e. Ci = in the linear alignment model). In the right hand panel the dotted (blue) and 
dot-dashed (magenta) lines show the inferred erg when the theoretical model varies 7 sca ic in only the GI or II terms. The gray 
band denotes the rough current uncertainty on erg from current data, of ±0.07. 



7scaie = 0.21. If we incorrectly analysed the data using 
the linear alignment model we would bias o% low by about 
0.1. This is relatively insensitive to the exact amplitude 
assumed for the linear alignment model, since the range 
of scales probed by the fiducial survey is dominated by 
the one-halo term (see Fig. ^ . 

An analysis that ignores intrinsic alignments corre- 
sponds to assuming a C± — in the linear alignment 
model, as shown by the leftmost point (cross) of the 
dashed (or dot-dashed) line in the left hand panel, for 
which er 8 is biased low by just under 0.1. This is because 
the fiducial ( "observed" ) model is contaminated by in- 
trinsic alignments, which are dominated by GI. There- 
fore the fiducial model power spectrum is lower than it 
would have been in the absence of intrinsic alignments. 
Attempting to fit to these "observed" points neglecting 
intrinsic alignments will result in an underestimate of cr g 
because the lensing-only power spectrum with the true 
a 8 will be too high. 

If the NLA model is used when fitting to the data, and 
yet if the true universe (fiducial model) followed the halo 
model, then the bias on as is reduced because the fitted 
model contains some of the suppression due to the GI 
effect. The bias is reduced significantly as the amplitude 
of the intrinsic alignment model is increased because the 
NLA mimics to a large extent the shape of the halo model 
on small scales. Fig. [9] shows that the NLA amplitude is 
generally smaller than that of the halo model, therefore a 
larger than fiducial C\ is required for the NLA to match 
the halo model. 



On fitting the halo model to the fiducial halo model, 
the bias on as is relatively insensitive to the exact am- 
plitude of the two-halo term. This is for the same reason 
that the linear alignment model fit is shallow, because 
the two-halo term has a small weight for this survey. 

We have investigated the effect of the assumed am- 
plitude of the one- halo term on biases on crg. In the 
right hand panel of Fig. [Til] we show the impact of as- 
suming the wrong one-halo amplitude, as parameterised 
by 7scaie , the effective reduction in stick length of the 
single halo sticks. The full result (solid line) is relatively 
complicated so we also show the contributions to this by 
varying 7 sca ie in only the GI term (dotted line) or in the II 
term (dot-dashed line). The dependence of the one- halo 
GI angular power spectrum amplitude on 7 sca i is simply 
linear, and the bias on ag is roughly linear in the intrin- 
sic alignment perturbation, for small intrinsic alignment 
contributions. The dashed line shows that as the ampli- 
tude of the GI term is increased in the fitted model then 
the fitted as increases due to the over suppression of the 
fitted model. 

The dependence of the one-halo II angular power spec- 
trum on 7scaicis quadratic (C/ 7 oc7 sca i c 2 ) therefore the 
bias on as is also quadratic in 7 sca i G . The direction of the 
bias is opposite to the GI term alone, because fitting a 
model with too much II contribution will wrongly boost 
the total angular power spectrum, and as will have to be 
reduced in the fit to compensate. 

Taking the two effects together explains the full result 
shown by the solid line in the right hand panel of Fig. QJjl 
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in which 7 sca i e is varied self-consistcntly in both the GI 
and II terms. The net effect is that an underestimate for 
7scaie tends to cause an underestimate of ag. 

Finally we consider the case where the fiducial model is 
the NLA, and the fitted model is the halo model (dashed 
line). The bias is the same shape as that when the fiducial 
model is the halo model, but there is an offset (of the 
same size as the dashed line in the left hand panel at 
Ci = 5 x 1CT 14 (h 2 M Q Mpc 3 )" 1 ). We see that the bias 
is removed on using around half the effective stick length 
of the fiducial halo model. 

In both panels of Fig. [10] we compare the biases to the 
one-sigma Ug error bar of 0.07 obtained by Ref. [4{|. The 
largest biases are comparable to this error bar. However, 
we note that the range of scales used in the actual cos- 
mological parameter analysis often smaller than the full 
range used in the angular power spectrum Fisher analysis 
we use here, and therefore the actual biases from intrinsic 
alignments will tend to be smaller. 

V. CONCLUSIONS 

We have presented a new model for intrinsic alignments 
of galaxies. This is inspired by the halo model for the 
matter and galaxy distribution in the universe and by 
simulations that suggest that satellite galaxies point to- 
wards the centers of halos. We assume the universe is 
entirely made up of halos that contain (i) a single cen- 
tral galaxy with an orientation determined by the cur- 
vature of the large scale potential (the linear alignment 
model) and (ii) satellite galaxies distributed spherically 
according to the halo dark matter profile that are ori- 
ented pointing at the center of the halo. 

We have described several terms that contribute to the 
intrinsic-intrinsic (II) and lensing-intrinsic (GI) source 
power spectra. The two-halo term describes contribu- 
tions to the power spectra from pairs of galaxies in two 
different halos and can in general contain contributions 
from central-satellite correlations, satellite-satellite cor- 
relations and central-central correlations. We find that 
the dominant contribution comes from the central-central 
correlation, and this corresponds to the existing well 
known linear-alignment model. 

The one-halo term describes correlations arising from 
pairs of galaxies situated in the same halo. In our sim- 
ple model the central-satellite term is zero and the one- 
halo term is equal to the satellite-satellite contribution. 
This is the major result of this paper. The correlations 
between orientations of galaxies within the same parent 
halo cause a boost in intrinsic alignment power spectra 
on small scales. Small scales are not expected to be at 
all well described by the linear alignment model. 

We have considered the impact of random perturba- 
tions to the satellite galaxy alignments and found that 
it reduces the GI power spectrum by around a factor of 
0.2 and the II power spectrum by a factor of 0.2 2 . We 
have demonstrated an analytical method for calculating 



these power spectra based on a multipole expansion. We 
have also provided fitting formulae for the one-halo term 
containing 9 parameters for each of the II and GI power 
spectra. 

The intrinsic alignment contributions to cosmic shear 
angular power spectra are illustrated for a deep tomo- 
graphic survey and a medium-deep survey using a single 
redshift bin. As expected the biggest contributions come 
from the lowest redshifts and for the non-tomographic 
survey the largest effect is a suppression of power due to 
the GI alignments. 

We made a rough estimate of the possible impact of in- 
correctly using the wrong intrinsic alignment model for 
cosmological parameter constraints from a medium-deep 
survey We find that, if intrinsic alignments are ignored 
completely in the cosmological parameter analysis, and 
all available scales are included in the analysis, then es- 
timates of the matter power spectrum normalisation as 
can be biased low by approximately the current 1-er un- 
certainty. However, we note that cosmological parameter 
constraints are usually performed using a limited range 
of scales, and therefore will be much less sensitive to the 
small scale effects introduced in this paper and <r% will be 
much less biased. Investigation of how the bias depends 
on the exact analysis is beyond the scope of this work. 

Our implementation of the halo model we present is 
subject to a number of limitations. We assume all halos 
are spherical and all satellites point at the halo center, 
subject to a random misalignment angle. Future work 
could investigate the impact of using elliptical halos in 
which satellites are anisotropically distributed and the 
satellite orientations depend on their position within the 
halo [llj]. We assume there is only a single population 
of satellite galaxies, and a single population of central 
galaxies. However a more sophisticated approach could 
use a mixture of spiral and elliptical galaxies with dif- 
ferent alignment models for each. We have also assumed 
that the ellipticity of the satellite galaxies is indepen- 
dent of radius and independent of the mass of the parent 
halo [2(| . Finally we have assumed that the central galax- 
ies are aligned with their host halo, whereas in fact there 
could be additional randomness [24[ . 

The halo model implementation could be improved 
considerably in the future as more fitting functions be- 
come available from simulators. Due to this separation 
into one and two-halo terms, significant progress can be 
made using high resolution simulations of single parent 
halos instead of full-sized cosmological simulations. The 
fitting formula for the distribution of galaxy alignment 
angles from Ref. [lj| was very useful. It would be help- 
ful to in addition implement fitting formulae for (i) how 
satellite galaxy ellipticity depends on radius within the 
parent halo, and on parent halo mass and (ii) how the 
distribution in alignment angles changes with radius and 
parent halo mass. For future extentions to this halo 
model it would also be useful to have (iii) the distri- 
bution in alignment angles between galaxies and dark 
matter halos (when using N-body simulation results to 
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construct a model for galaxies) , (iv) fitting formualae for and 
how the satellite galaxy orientation depends on position 
within the halo for non-spherical halos, (v) the alignment 
distribution between satellites and the central galaxy in 
a halo, and (vi) a dependence of all the above on satellite 
mass and/or type. Once point (v) is known, a central- 
satellite one-halo term could potentially be added to our 
results using the distribution between the ellipticity of 
the central region of a halo and that of the total halo 
recently given in Ref. [24| . To build a halo model for f j^ e f 
central galaxies that could replace the LA model used 
here for the two-halo term, one would also need to know 
the IA correlation function between elliptical halos. We 
will pursue this in future work. 
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Again following Ref. [51|, if we assume u>{9) <C 1, 

co(0)/(l + ui(9)) ui(9). If we further approximate oj(9) 

as a power law in 6, then the integral in Eqn. (|A6[) can 

be done with Eqns. 11.4.33, 11.4.34 in Ref. [II]. We use 

the (sli ghtl y modified) power-law model from Eqn. (66) 
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APPENDIX A: NORMALIZATION WHEN 
SOURCES ARE CLUSTERED 

Because galaxy positions are correlated, we must nor- 
malize the ellipticity power spectra from Section IIIII by 
the number of galaxy pairs at a given angular separation, 

N paila (9) = n 2 ^(l + u(9)), 

where u)(0) is the galaxy angular correlation function. 
Applying this normalization to the ellipticity correlation 
function, we get, 



e ± (*,x)=/ ^ JoA 



o 2?r 1 + u(9) 



p ee (£,x)±p b i b (£,x) 



(Al) 

Following a similar derivation in Ref. [5 II ] - we can rear- 
range the multiplicative normalization into an additive 
correction using, 



1 



g(g) 



(A2) 



To get the additive correction to the ellipticity power 
spectra, we put Eqn. (|Aip into the transformation, 



p EE ' BB (i) = nj ede [£ + (0)Mm ± Z-{e)Ute)] . 

(A3) 

to get the modification to the power spectra due to the 
extra u(6)/(l + w(6), 



with 7 = 1.7 and A = 0.5. 

With this model, we find negligible corrections to the 
unnormalizcd spectra given in Section IIIII and thus we 
neglect the effect of u(9) in Eqn. (|A1[) . 



APPENDIX B: COMPUTATION OF 
ELLIPTICITY POWER SPECTRA 



In this section we first derive the exact dependence of 
the satellite ellipticity model on the radial alignment an- 
gles (3 and r\ from Eqns. (fH)|) and (TT|) and then show how 
we compute the power spectra by means of a multipole 
expansion of the satellite cllitpicity distribution in a halo. 



1. Galaxy major axis as a function of radial 
alignment angles 



First consider the explicit dependence of the ellipticity 
of the satellite galaxies on the alignment angles (3 and r\ 
from Section fill Bl If we write 

e = (sin# e cos <p e , sin(9 e sin </> e , cos 9 e ) , (Bl) 

then the ellipticity of a satellite galaxy in a halo is, 



1 r°° 

ap e i ebb (£, x ) = -t: / e'd£'{j I0 (e,i')[p EE (£') 

7 2 J Q 

+P BB (£')] ± J M {£,£') [P EE {£') - P BB {£')}}, (A4) 



with 



AP SAl (£, x )^2n e'de'P StV (£', X )J I2 (£,i'), (A5) 
Jo 



j(t) = j(r, m) sin0 e (cos 20 e , sin 2^> e ) . (B2) 

To find 9 e and <p e as functions of 9, </>, f3, and rj, we use 
a set of Euler angles to rotate to a coordinate system r' 
with z 1 = f. This can be accomplished by setting the 



Euler angles (in the zxz convention) to a E 
f3 E = }J E = -tt/2. 



+ 7T/2, 



Explicitly, the components of the unit vector describing the major axis of a galaxy in the 2 coordinate systems are 
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related by 



e = (sin (3 cos 77, sin (3 sin 77, cos (3) T 



R(a E 7 f3 E }1 E )e 




1 

cos f3 E -sin/3^ | [ sina B cosa B | e. (B3) 
sm/3 B 




We then have e = i? x e' and 



sme/ e = w ef + e 2 ,, cos< 



ei 



'■2 



(B4) 

giving sin# e , cos</> e , and sin</> e as functions of 0, (f>, /3, 
and 77 as desired. 



2. Satellite ellipticity density run 

We define the normalized density run of the satellite 
ellipticity distribution in a halo as 



zk-r 



w(k\m 



J d 3 rj I (r,m) p NFW (r, m,c) e 
/ d 3 r |7 / (r, m)| p NFW (r,m, c) 



(B5) 



where, when allowing for a distribution in radial align- 
ments, 

I 1 (r, m) = 7(777) sin (0 e (0, 0, (3, 77)) exp [2i0 e (0, <j>, /3, 77)] , 

(B6) 

where e and e describe the orientation of the galaxy 
major axis as in Eqn. (|B1[) . and we have assumed that 
the magnitude of the intrinsic ellipticity, 7(771), is inde- 
pendent of position within the halo. Separating the two 
components of the complex ellipticity, 



7«i(k|m) = — J d 3 r cos(20 e ) sm(d e )u(r\m)e lk r , 

w 2 (k\m) = ^- J d 3 r sin(20 e ) sin(0 e )u(r|TO)e jk - r . (B7) 

We can simplify the computation of this transform by 
performing a multipole expansion, 



;(k|m) = J2(2£ + 1) i e u e (k\m) 
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OO 
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where cos 7 = sin0fc sin (9 cos( 



cos 61. cos0 and 



Ug(k\m) = / r dr u(r\m) je(kr) 
Jo 



are normalized multipole moments of the NFW profile. 
The integral over <j> in ft guarantees that ft = Q for odd 
£ (and I = 0). The first four even multipolcs ue(k\m) 
are shown in the left hand panel of FigHT] If we as- 
sume that the radial mis-alignment of the galaxies is 
uncorrelated such that for 2 galaxies P(/3i, 771, /3a, 772) = 
P(/3i_, 771) P(/32, »72), then this expression for /g can be in- 
tegrated over /3 and 77 at this stage before being used to 
compute correlation functions. 

For the case of j3 = (so 8 e — 9 and <fi e = 4>), all the 
integrals over <fi in are of the form, 



cos 2(j> 
sin 2(f) 

with, for example, 



cos 



.92 = .94 



96 



157T 

~32~ : 



.98 



7tt 
16' 



5io 



cos 2cj) k 
sin 2(j) k 



1057T 



256 



9n 

(B9) 
(BIO) 



This makes it clear that the only (f> k dependence in ft is 
in the cos/sin phase factors. 

We can then use Eqn. @ to write, 

7 £(k) = (cos 2 (20 fc ) +sin 2 (2^)) ^(k)! 
= |7 J (k)| 

7&(k) = (cos(20 fc )sin(2</> fc ) -cos(20 fe )sin(20 fe )) | 7 7 (k)| 
= 0. (Bll) 

So, the stick model with perfect radial alignments for 
satellite galaxies has zero B-mode. 

Again for /3 = 0, the angular term in the multipole 
expansion, ft, can be computed using the formula, 



ll\Pki 



= e Pi 

m=0 

x sin J (0fc) cos 



E 

J'=0 



.7 



.9/ 



777 — J j 

(B12) 



where p; m are coefficients of the £th Legendre polynomial 
and 



/(a, 6) 



(1 -x 2 ) 2 x b . 



(B13) 



We compare Eqn. (|B8|) for £ < 10 with the exact com- 
putation of Eqn. (|B5[) in the right hand panel of Fig. QTJ 
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k [h/Mpc] 



FIG. 11: The left hand panel shows the first 4 nonzero multipole moments of the normalized NFW profile. Higher multipoles 
are subdominant for small wavenumber. The right hand panel shows a comparison of the Fourier transform of 7 7 for the stick 
model for satellite galaxies with perfect radial alignments computed using a 3-D FFT (red, solid) and the £ = 2, 4, 6, 8, 10 
multipoles of the NFW profile. 



We find that the sum over multipoles converges rapidly 
on the scales of interest. All of the results in the main 
body of the paper use only the £ = 2 term when comput- 
ing ui(k|m). 
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